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Abstract — This work will investigate the behaviors of Moving distributed masses of orthotropic rectangular 
plates resting on bi-parametric elastic foundation with clamped-clamped end conditions. The governing 
equation is a fourth order partial differential equation with variable and singular co-efficients. The solutions to 
the problem will be obtained by transforming the partial differential equation for the problem to a set of coupled 
second order ordinary differential equations using the technique of Shadnam el al[18] which shall then be 
simplified using modified asymptotic method of Struble. The closed form solution will be analyzed, resonance 
condition shall be obtained and the result will be presented in plotted curves for both cases of moving distributed 
mass and moving distributed force. 

Keywords — Bi-parametric foundation,orthotropic, foundation modulus, critical speed, resonance, modified 
frequency. 


I. INTRODUCTION 

A plate is a flat structural element for which the thickness is 
small compared with the surface dimensions or a plate is a 
structural element which is thin and flat. By "thin", it means 
that the plate’s transverse dimension, or thickness, is small in 
comparison with the length and width dimensions. That is, the 
plate thickness is small compared to the other dimensions. A 
mathematical expression of this idea is: 

^«1 ( 1 . 1 ) 

where 'T' represents the plate’s thickness, and ’L’ 
represents the length or width dimension. The thickness of a 
plate is usually constant but may be variable and is measured 
normal to the middle surface of the plate. Plates subjected to 
in-plane loading can be solved mostly by using 
two-dimensional plane stress theory. On the other hand,plate 
theory is concerned mainly with lateral loading. There are 
whole lot of differences between plane stress and plate theory. 
One of the differences is that in the plate theory,the stress 


components are allowed to vary through the thickness of the 
plate,so that there would exist bending moments. 



Fig.1.1: The Skeletal Structure of Orthotropic Plate 


The vast majority of the world’s bridges with record long 
spans have utilized orthotropic steel deck systems as their 
superstructures.These types of decks have been used 
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extensively in Europe, Asia and South America. 
Approximately 6,000 orthotropic steel deck bridges exist in a 
world of about six million bridges. With the growing trend 
quicker construction practices with an overall longer bridge 
life. The other leading benefits of this bridge decking system 
are the minimization of dead load in the design and the rapid 
construction that will lessen the impact on traffic. 

Tracing the plate theory to its roots, one travels back to the 
American revolution time. During this revolutionary period, 
several scientists and engineers performed numerous 
researches on plates. Euler [1] performed a free vibration 
analysis of plate problems and indicated the first impetus to a 
mathematical statement of plate problems. Chladni [2], a 
German physicist, performed experiments on horizontal 
plates to quantify their vibratory modes. He spattered sand on 
the plates, struck them with a hammer, and observed there 
were regular patterns formed along the nodal lines. Bernoulli 
[3] then attempted to theoretically justify the experimental 
results of Chladni using the previously developed 
Euler-Bernoulli bending beam theory, but his results were 
unable to capture the full dynamics. The Erench 
mathematician Germain [4] developed a plate differential 
equation that lacked a twisting term but one of the reviewers 
of her works, Lagrange [5], corrected Germain’s results. 
Thus, he was the first person to present the general plate 
equation as stated by Ventsel and Krauthemmer [6]. Cauchy 
[7] and Poisson [8] developed the problem of plate bending 
using general theory of elasticity. Poison successfully 
expanded, the Germain-Lagrange plate equation to the 
solution of a plate uder static loading. In this solution, the 
plate flexural rigidity was set to be a constant term D, as 
showed in Ventsel and Krauthammer [6]. Navier [9] 
considered the plate thickness in the general plate equation as 
a function of rigidity D. Kirchhoff [ 10] shed more light on thin 
plate theory based on some assumptions that are now referred 
to as "Kirchhoffs hypothesis". These hypotheses are 
fundamental assumptions in the development of hnear, 
elastic, small-deflection theory for the bending of thin plates. 
It is based on Kirchhoff hypothesis that straight lines normal 
to the undeformed mid-plane remain straight and normal to 
the deformed mid-plane. In accordance with the kinematic 
assumptions made in the classical plate theory, all the 
transverse shear and transverse normal strains are zero. In the 
works of these great scientists, the flexural rigidity was 
considered to be constant, D. That is to say, they all treated 
plates as homogeneous and isotropic materials, that their 
material properties remain unchanged in all directions but in 


reality plates are orthotropic. 

In the 1920’s, American engineers began using steel plate 
riveted to steel beams for large movable bridges. The purpose 
was to minimize the dead load of the lift span. In 1938, the 
American Institute of Steel Construction (AISC) began 
publishing reports on the steel-deck system. AISC called this 
the "battledeck floor" because it felt the steel deck had the 
strength of a battleship. The orthotropic deck was a result of 
the "battledeck floor". This floor consisted of a steel deck 
plate, supported by longitudinal (normally I-beam) stringers. 
In their ton, these stringers were supported by cross beams. 
Eollowing World War II, German engineers developed the 
modern orthotropic bridge design as a response to material 
shortages during the post-war period. The orthotropic deck 
reduced the weight of continuous beams considerably and 
permitted spans and slenderness ratios unknown until then. 
Since the 1950s, steel highway bridges have been constructed 
with a steel deck structoe. This structure is called orthotropic 
deck. These decks were constructed with a deck plate 
supported by stiffeners of various shapes and by cross-beams 
and main girders. The deck plate acts as the top flange of the 
deck girders. 

Since 1965, a new generation of orthotropic decks has been 
used, with cold-formed trapezoidal stringers, so-called 
troughs. This trough permitted cross-beam spacing. An 
orthotropic deck or orthotropic bridge is one whose deck 
typically comprises a structural steel deck plate stiffened 
either longitudinally or transversely, or in both directions. 
This allows the deck both to directly bear vehicular loads and 
to contribute to the bridge structure’s overall load-bearing 
behaviour. The orthotropic deck may be integral with or 
supported on a grid of deck framing members such as floor 
beams and girders. Decks with different stiffnesses in 
longitudinal and transverse directions are called ’orthotropic’. 
If the stiffnesses are similar in the two directions, then the 
deck is called ’isotropic’. The steel deck-plate-and-ribs 
system may be idealized for analytical piuposes as an 
orthogonal-anisotropic plate, hence the abbreviated 
designation ’orthotropic.’ In vibration analysis of plates, the 
effect of elastic foundation must be considered when plates 
are mounted on elastic springs such as pavement of roads, 
footing of buildings and bases of machines.Extensive 
understanding of vibration behavior of orthotropic plates is 
important in many engineering fields. Therefore, the vibration 
analysis of plates resting on elastic foundation has been the 
subject of many researches. Ugural [11] and Okafor and 
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Oguaghamba [12] applied Ritz and Galerkin’s methods to 
solve isotropic and orthotropic plate problems. Their results 
are very close to those exact solutions. Some researchers 
applied superposition method for analysis of vibrations of 
plates having classical boundary conditions for instance Ohya 
et al [13] for deformable plates.Awodola and Omolofe [14] 
investigated the response of concentrated moving masses of 
elastically supported rectangular plates resting on Winlker 
elastic foundation by the method of variable separable. 
Agarana, Gbadeyan and Ajayi [15] solved an isotropic elastic 
damped rectangular Mindlin plate resting on Pasternak 
foundation by using numerical method. Are, Mown and 
Gbadeyan [16] dealt with a damped simply supported 
orthotropic rectangular plates resting on elastic Winkler 
foundation by applying method of separation of variables to 
reduce the problem to second order coupled differential 
equation and then applied numerical methods to obtain the 
solution. Gbadeyan and Dada [17] improved on their previous 
works by considering the dynamic response of a Mindlin 
elastic rectangular plate subjected to distributed moving load, 
but neglected the effect of damping. 

The orthotropic materials are of a great importance and 
interest in fields of the modern industrial technologies 
applications due to their hardness,lightness and super 
elasticity.The problem shall be solved analytically, which is 


not found in the work of the early researchers by using the 
technique of Shadnam et al[18]. Barlie [19] carried out the 
thermo-mechanical morphological study of CFRP in different 
environmental conditions. Guo [20, 21] examined the fatigue 
lif of reinforced concrete beams strengthened with CFRP 
under beading load.Hosseini [22] investigated the roadway 
metallic bridge using non-prestressed bonded and prestressed 
unbonded CFRP plate. In the same view, Ju [23] and Li [24] 
examined the fatigue life of reinforced and non-prestressed 
bridge desk under variable amplitude load and different 
adhesives. Liang [25] studied the reliability analysis of bond 
behavior of CFRP concrete under wet-dry cycles. Xin Yuan 
[2] also investigated the fatigue performance and life 
prediction of CFRP plate. 

II. GOVERNING EQUATION 

The dynamic transverse displacement V(x,y,t) of 
orthotropic rectangular plates when it is resting on a 
bi-parametric elastic foundation and traversed by distributed 
mass Mj. moving with constant velocity along a straight 
line parallel to the x-axis issuing from point y=s on the y-axis 
with flexural rigidities and Dy is governed by the fourth 
order partial differential equation given as 


D. ^ W[x. y. t) + 2B Vl/(x, y, t)-I-Dy^ y, t) + /r ^Wix, y, t) - phR^ 

^ gjgf .2 ^V’ + gylgt 2 ^ y’ y. 0 - ^0 (x, y,t)-\-^W(x, y, t)] 

- ^ [MrPHix - ct)H(y -s)-M^(^W(x, y, t) + 2c^ W(x, y, t) + c?-^W(x,y, t) 

)H(x - c^t)H(y - s)W(x,y, t)] = 0 (2.1) 


where Dj. and Dy are the flexural rigidities of the plate along x and y axes respectively. 


Ex and Ey are the Young’s moduli along x and y axes respectively, Gg is the rigidity modulus, and Vy are Poisson’s ratios 
for the material such that E^Vy = EyV^ , p is the mass density per unit volume of the plate, h is the plate thickness, t is the time, x 
and y are the spatial coordinates in x and y directions respectively, Rg is the rotatory inertia correction factor. Kg is the foundation 
constant and g is the acceleration due to gravity, H(.) is the Heaviside function. 


Rewriting equation (2.1), one obtains 
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li^.Wix,y, t) + fiCP^Wix,y, t) = phR,Wix,y, t) + Wix,y, t)] - 2B 

Wix, y, Wix,y, t)-Dy£-, Wix,y, t) - K,Wix, y, t) + pco^Wix,y, t) 

+ Go[^Wix,y,t) + ^^Wix,y,t)] + T,?=i [M^gHix - c^miy - s) - M,(^M/(z,y,t) 

-\-2c^^Wix,y,t) + ^wix,y,t))Hix - Cr-tjHiy - s)Wix,y,t)] (2.2) 

which can be expressed further as 

^Wix,y,t) + colWix,y,t) = t) + ^Vl/(z,y,t)] - f ^ 

Wix,y, t) - ^^Wix,y, t) - Vl/(z,y, t) + K - jWix, y, t) + ^ 

VK(z,y,t) + ^Vl/(x,y,t)]+2:!y=i [^Hix - c,t)Hiy - s)-^i^,Wix,y,t)2c 
^Wix,y,t) + cj^,Wix,y,t))Hix - c^miy - s)Wix,y,t)] 


where co^ is the natural frequencies, n = 1,2,3,... The initial conditions, without any loss of generality, is taken as 
Vl/(x,y,t) = 0 = ^Vl/(x,y,t) (2.4) 

III. ANALYTICAL APPROXIMATE SOLUTION 

In order to solve equation (2.3), one applies technique of Shadnam et al[19] which requires that the deflection of the plates be ii 
series form as 

Wix,y,t) = 0„(z,y)(3„(t) (3.1) 

where 0„(:r,y) = 0„i(z)0„j(y) and 

0„j(x) = sin-^z + Ajiicos-^x + B„jSinh-^z + Cnicosh-^x 


^ni(y) = sin^y + A„yCOS^y + B„jSinh^y + C„jCOsh^y (3.2) 

The right hand side of equation (2.3) written in the form of series takes the form 

2:n=i /?o[^Vl/(x,y,t) + ^Vl/(x,y,t)]-^^Vl/(x,y,t)-^^Vl/(x,y,t) 

+ ^^n-'jW(x,y,t) -{-^[^Wix,y,t) -{-£^Wix,y,t)] + 

[^^riM^iy - s) -^i^,Wix,y,t) + 2c,£-Wix,y,t) + £ 

Wix,y,t')')Hix - CytjHiy - s')Wix,y,t')] = ‘I’n(^-y)o'n(f) 

On multiplying both sides of equation (3.3) by d5m(x,y), integrating on area A of the plate and considering the orthogonality of 
Om(^-y)> one obtains 
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^nit) = L - 

7 £ M/(x, y, t) - ^^ M/(x, y, t + (col - ^)Vl/(x, y, 0 + ^ W^x, y, t) 

+ M/(x,y, t)) + - Crt)H(y - 5 ) - ^ Vl/(z,y, t) + 2c, ^ 

W(x,y,t') + c^ ^ VK(x,y, t')')H(x — Crt')H(y — s')W(x,y, 


and zero when n ^ m 
where 

e* = <^l{x.y)dA (3.5) 

Making use of equation (3.1) and taking into account equations (2.3) and (3.4), equation (3.3) can he written as 

0„(x,y)[Q„(t) + (plQnm = 2^2]“ 1 [/?o(^^^0^(x,y)(3,(t) + 

0„^(x,y)(3,(t)) -y^j5S^O^(x,y)(?,(t) -^^^^0„(x,y)(3,(t) 

+ {col - ^)0,(x,y)0^(z,y)(3,(t) + ^(^^^^0„(x,y) 




0„^(x,y)Q,(t) + 2c,^^|^0^(x,y)(3,(t) + c,2^^^^0^(x,y)(3,(t)) 
//(x - Crt)H{y — s))]d4 


(y-5) -^(%(x,y) 
(3.6) 


On further simplification of equation (3.6), one obtains 

_ 1^ r 5^0„(x,y) .. d^O„(x,y) 

Quit) + iolQnit) = - 2 ^ J [fio( g ^2 ^mix,y)Q,it) + —cp^(x,y)Q,(t) 


2g d"Oq(x,y) 
/r 5^2 5y2 


D,,d*^Jx,y) Dyd^^Jx,y) 

.ix,y)Q,it) - -j 0^(x,y)Q,(t) --j 0^(x,y) 


, ffo Go 520„(z,y) 52(ho(x,y) 

0,(0 + - y)‘i’,(^<y)‘i’m(^<y)0,(0 + y( <i^m(^.y)0,(0 + 

‘I’m(^<y)0,(O) + y (^ 0 „(x,y)//(x - c,t)//(y - 0 - ^( 0 ,(x,y) 0 „(x,y)( 3 ,(t) 

^ fi P 


+2c/-^^^^mix,y)Qcit) + c^^-^^^mix,y)Qcit))H(x - Crt)H(y - s))]dA 


The system in equation (3.7) is a set of coupled ordinary differential equations 

Using the Fourier series representation, the Heaviside functions take the form 


H(x - c^t) =2 + -^r 


sin(2n+l)7r(3:-c,t) . 


H(y-s) = ^-+, 0 < y < 1 

On putting equations (3.8) and (3.9) into equation (3.7) and simplifying one obtains 


http://www.aipublications.com/iicmp/ 


Page I 75 
















Open Access 


International journal of Chemistry, Mathematics and Physics (IJCMP) 
[Vol-4, Issue-4, Jul-Aug, 2020] 
https://dx.doi.ora/10.22161/iicmp.4.4.2 

ISSN: 2456-866X 


Quit) + co^Qnit) - [RoToQ.it) - jT,Q,it) - 4 

(col - + ^fT,Q,(t) - ^ ((Te + ^ (XT=i El 

, El , ^3* El ^ 


^ 2;+l 

^voo cos(2; + l)7rc, 

C^y=i *5 ^77^“ 


"‘ = ^ 2fcTI--^^=1 - 2fc+l 

_yoo sin(2j+l)7rCyt >. 1 ^yoo cos(2fe+l)7r5 

2y+l ^ 2y+l ^ 2fc+l 

1 cosr2/+l')7rc^t 


2fc+l 

_ ^^J-))Q,{t) + 2Crt(T, 

sin(2;+l)7rc. 


,,,, ....—1:7=1^10 

'rt^.voo I-* cos(2fc+l)7rs „oo i-* sin(2fc+l)7rs- 1 ,„oo i-* 

.... -)Sfc=l ^11 -2fe=l ^12 2fc+l ) + = 1 ^13 

cos(2; + l)7rc.t _ sin(2; + l)7rc.t J_ cos(2fe+l)7rs _ 

2j + l 2; + l 2fc+l ^k-1 16 

sin(2fc+l)7i:S '' . (f^ + r^cp i i_ .y<?> g* cos(2; + l)7rCrt „oo sin(2; + l)7rCrt 

2fc+l r'. 8 JJ.2 vZ,j = l 17 2; + l ^J = l 18 2;+l 

^.Voo c* cos(2k+l)7rs „oo r.. sin(2k+l)7rs- 1 ,„oo r-* cos(2;+l)7rCrt 

)ak=l il 9 Lk^l E 20 2 k+l ^ + to <'^ 7=1 ^21 

„oo I-* sin(2; + l)7rCrt. 1 ,„oo i-* cos(2k+l)7rs „oo i-* sin(2k+l)7rs-. 

- Zj=l E22 --) + - (Zk=l E23 -2:^=1 ^24 2k+l )) 


which is the transformed equation governing the problem of an orthotropic rectangular plate resting on bi-parametric elastic 
foundation. 


To = [il^%(.^>y)^mix,y)-^^%(x,y)<t>^(x,y)]dA 

(3.11) 

Ti = ^[^M^,y)]^mix,y)dA 

(3.12) 

Ti = I,i^[^q(x,y)]^m(x,y)dA 

(3.13) 

T3 = I,,£-A<i>q(^'y)]'i>mix,y)dA 

(3.14) 

Ta = ^q(x,y)^m(X,y)dA 

(3.15) 

Ts = [-^^qix.y) + -^^q(x,y)]^m(x,y)dA 

(3.16) 

T6=reh'^<i^^’y^^rn(x,y)dA 

(3.17) 

El = <P^(x,y)<Pm(x,y)sin(2j 4- l)7rzd.4 

(3.18) 

El = '^q(x,y)<t>m(x,y)cos(2j 4- l)7rzd.4 

(3.19) 

El = Oq(x,y)Om(x,y)sin(2fe 4- l)7ryd.d 

(3.20) 

El = Oq(x,y)Om(z,y)cos(2fe 4- l)7ryd.d 

(3.21) 

II 

II 

(3.22) 

Tj = il^q(X’y)^m(X,y)dA 

(3.23) 

El = ^(0q(x,y))0m(x,y)sin(2j 4- l)7rxd.4 

(3.24) 

^10 = ^ ('^qCX' y))'^m(x, y)cos(2j 4- l)7rxd.4 

(3.25) 
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£■ 1*1 = ^(‘I’q(^-y))‘I’m(^-y)sin( 2 fc + l) 7 rydy 4 

(3.26) 

^12 = ^‘I’q(^-y)‘I’m(^-y)cos( 2 fc + l)7rydi4 

(3.27) 

El-i = Eg, Bi‘4 = Eio, Eis = Eji, Elf, = EI 2 

(3.28) 

^8 = ^i^qi^’y))^m{X.y)dA 

(3.29) 

Ell = ^i^qix,y))^rn(x,y)sm(2j + l)7rxdA 

(3.30) 

E*g = ^(Oq(x,y))Om(:r,y)cos( 2 ; + l)7rzd.4 

(3.31) 

Big = Oq(z,y)Om(^.y)sin( 2 fc + l)7rydy4 

(3.32) 

£■20 = ^(‘I’q(^-y))‘I’m(^.y)cos( 2 fc + l)7ryd^ 

(3.33) 

E 2 I — Ely, E 22 = Big, B 23 = Big, B 24 = B 20 

(3.34) 


is assumed to be the products of functions ^pm(^)^bm(.y) which are the beam functions in the directions of x and y 
axes respectively. That is 


= (3-35) 

where 


‘I’bm(y) = sinA„m + AtmCOsAprny + BbrnSinhA^my + CtmCOShAp^y (3.36) 


where Bp^, Cp^, Ap^, and are constants determined by the boundary conditions. And Op^ and are 
called the mode frequencies 
where 


(3.37) 


Considering a unit mass, equation (3.10) can be rewritten as 


Quit) + colQnit) - 1 [RoToUt) - 

+ (0)^ - + ^(S” 1 Ej - 


yioo p+ sin(2y+l)7rct '. ,y,oo cos(2fc+ l)7rs 

■^f=l 2 2j + l AZ,fe=l 3 2fc+l 

cos(2; + l)7rct _ „oo p* sin(2; + l)7rct , J_ .v 

2;+i 2^; = 1^6 2;+1 

sin(2fc+l).s . + i-rV” , 




i:7=r 


i (j;- £cos(2^^ _ 

2;-H ^ IS 2fc+l ^k-1 16 


cos(2j + l)7rct 


2J+1 

sin(2k+l)7rs 

2k+l 


))Q,(t) 
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-2:”= i £2 o - 


^( 2 : 7=1 ^21 


cos(2j + l)7rct _ 


2 : 7=1 ^22 


+ ^(E7=i 2 : 7=1 

= 2 : 7=1 ^‘I’m(Ct)0^(5) 

equation (3.38) is the fundamental equation of the problem, where 

ur = ^, <p=L^L, 

O^(ct) = sina^(t) +i4^cosa^(t) + B^sinha^(t) + C^cosha^(t) 
O^(5) = sinl„ + .4^cosl, 


. =£=t£, 1 =E=t£ 

™ Ly 

3.1 Orthotropic Rectangular Plate Traversed by a Moving Force 

In moving force, we account for only the load being transferred to the structure. In tl 
0 = 0 in the fundamental equation (3.38), one obtains 

Quit) + (1 - -j^[pRoToQn(t) - 2BTMt) - DJMt) - DyTMt) 

-KoT.Qnit) + GoTsQnit) + E7=i,,^n ipRoToQcit) - 2BT,Q,it) - DJ^Q,{t) - Dy 
T,Q^(t) + iiicol - KJ,)Q^{t) + = 2‘I’m(ct)0„(s) (3.43) 

which can further be simplified as 

Quit) + - y[/tflo7’o<?n(t) - 2BTMt) - DJ^Qnit) - DyT^t) - KJMt) 

+GJ^Q^{t) + E7=i.,^n {pRJ^Qdit) - 2BT,Q^(t) - DJ^Q^it) - DyT,Q^(t) + 

- WQ,it) + GJsQ.im = yMgcp^ctjcp^is) (3.44) 

where = (1 - ^)o)^ 

Expanding and re-arranging equation (3.44), one obtains 

[1 - YBRoTo]Qnit) + (f^ - yJe)Qnit) - 72 : 7 = 1 , ipRoToQqit) - 2BT,Q^it) - 
Dj2Q,{t) - DyT.Q^it) + {puyl - + GJ^Q^{t)) = yMg^^{cm^{s) 


(3.39) 

(3.40) 

(3.41) 

(3.42) 


the inertia effect is negligible. Setting 


Simphfying 


Quit) 


j!tyhL 


Qn(t) 4 


[i-rMfioro] 

Q,(t) - DyT^Qq(t) + (/im^ - KoT^)Qq(t) + GoT^Q^^t)) 


E” (gRoToQ,{t) - 2BT,Q^{t) - D,T, 

brT‘I’-(ct)0™(5) 


= Je = -2BT, - DJ 2 - DyT, - KJ, + GJ^ 


For any arbitrary ratio 7 , defined a: 


I' =r 
7* 

l-y* 

For only 0 ( 7 *), 


7 = 


7*-f o(y*2)-f... 
obtains 
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y = y* 


On application of binomial expansion, 

i-rWo " ^ + o(y‘")+- ■■ 

On putting equation (3.48) into equation (3.46), one obtains 

Quit) + (fn - y76)(i + y*pRoTo + o(y*7+-■ ■ )<?n(0 + y*(i + y*pRoTo + o(y*^) 

(pRoToQ,(t) - 2BT,Q^(t) - D,T,Q^(t) - DyT,Q^(t) + - (3.49) 

WQcit) + GoTsQ.it)) = y*aMgil + y*gR,To + o(y*^)+...)0^(ct)0^(5) 

Retaining only o(y*), equation (3.49) becomes 

Quit) + [f^(i + y*pRoTo) - y*Je]Qnit) + y* (pRoToQqit) - 2BT^Q^it) - 

- DyT^Q^it) + - KoT^-)Q,(t) + GoTsQ,(t)) = y*Mg^Uct)<t>m{s) 

which is simplified further as 

Quit) + + y‘ E”=i,,^n (MRoToQqit) - 2BT^Q^it) - - DyT^Q^jt) 

- KoT^-)Q,(t) + GoTsQqCO) = y*Mg<t>^(ct)^rn{s) 

where 

j7 = [aii+y*BRoTo)-yVe] 

Using Struble’s technique, one obtains 

fnn = fn-(^) 

which is the modified frequency for moving force problem. 

Using equation (3.53), the homogeneous part of equation (3.51) can be written as 
Quit) + ^LQnit) = 0 


(3.48) 


(3.50) 

(3.51) 

(3.52) 

(3.53) 

(3.54) 


Hence, the entire equation (3.51) gives 

Quit) + ^LQn(t) = y‘M50„(ct)0™(5) 


On solving equation (3.55), one obtains 


Quit) = 


?nnC<-?in) 


[(am + fnn)(“mSinf„„t - f„„sinamt) - -4mfnn(“m + fnn) 


(cosa^t - COSf„„t) - BjnioCm - fnn)(“mSinf„„t - f„„sinhamt) + Cmfnn(“m 
-fnn)(C0Shamt - COSf„„t)] 


which on inversion yields 


(3.55) 


(3.56) 


Wix,y, t) = Epm=l E”m=l - f„„sinamt) - 

fnn(«m + fnn)(C0Samt - COS^^t) - Rm(«m “ - 4„sinhamt) + 

^nnioCm “ fnn)(coshamt - cosf„„t)](sin^x + A^^cos^-^x + RpmSinh^z + Cpm (^.57) 

cosh ^ z) (sin + A^^cos^-f^y + B,mSinh^y + C,mCosh^y) 

which is the transverse displacement response to a moving force. 

3.2 Orthtropic Rectangular Plate Traversed by a Moving Mass 

In moving mass problem, the moving load is assumed rigid, and the weight and as well as inertia forces are transferred to the 
moving load. That is the inertia effect is not negligible. Thus m ^ Q and so it is required to solve the entire equation (3.38).To 
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solve the equation,one employs analytical approximate method. This method is known as an approximate analytical method of 
Struble. The homogeneous part of equation (3.38) shall be replaced by a free system operator defined by the modified frequency 
Thus, the entire equation becomes 


Quit) + ^LQn(t) + [(T, + 1 E, 

sin(2y+l)7rct ' /'yoo p* cos(2fc+1)7 




2j+l 2fe+l 

cos(2j+l)7rct ^00 r,* sin(2;+l)7rct> 

^7^1-- 

sin(2fe+l)i 
2k+l 


+ 2c{T, 


-1:^=1 El 

'-) + ^{Y,u=iE;^- 
: 2 ;+i) 

2; + l 


,, C0S(2; + l)7rct 
2J+1 

in(2fc+l)7rs 


-) + i 


2fe+l 

cos(2k+i)ns 
2k+l 

1 /'V’oo cos(2y+l)7rct ^00 

= l 1^9 2 J + 1 


0:7=1 El 

-1.^=1 El 

f,, sin(2j+l)7ret 
^10 2/+1 


)0:“=1 ^1*1 - 2:7=1 ^1*2 + ^0:;” 1 ^1*3 

voo 17+ sin(2j+l)nct. 1 ^y’oo c* cos(2fc+l)7rs y’oo c* sin(2fe+ 
“ = l 2; + l ) + ^ (^^=1 ^15 -ik=l ^16 

))Ut) + c\T, + ( 2 : 7 =, _ 2 -, E,*, 


--2:7=1 

))Q,(t)] = 2 : 7=1 

where V* — 

On expanding equation (3.58), one obtains 

Quit) + ^nQuit) + rncp*[{T, + ^ (2:7=1 Fi* —^.. 

,.,^00 cos(2k+l)7rs „* sin(2k+l)7rs, 1 ,„oo r.* 

)ak=i £3 -lk=i £4 2fc+i ) + i;;(2.1=1 


- 2 : 7=1 ^ 2 ^ 


,, cos(2; + l)n:ct 

27+1 

- 2:7=1 El ± (J:-^ e ; £2£gl^ _ 5;-^ 

)Q „(0 + 2 c(T, + ^ ( 2:711 e; - 2711 Fi*o 


(3.58) 
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^OO t 7* Sin(2; + l)7rct^ , 1 cos(2k+l)7rs r-* sin(2k+l)7rs 

- ^;=i ^14 iTTi ; + ^ ak=i iis Lk=i ii6 2k+i 

))0n(t) + + ;^ ( 2 : 7=1 ^1*7 _ 27 .^ F,*, 

(-voo r-* cos(2k+l)7rs sin(2k+l)7rs^ , 1 r-* cos(2j+l')KCt 

^20 2k+l ^ + to = l ^21 ^77^^ 

v,oo n* sin(2; + l)7ret^ , 1 cos(2fc+l)7rs sin(2fc+l)7rs 

- Lj=l E22 2 ; + l ) + to ^23 ^777 -lk=l ^24 2 fe+i 

))(3„(t)] [(7, +^(2”, 

Sin (2;+l)7rct ^^^^ cos(2fc+l)7rs . 


n- - ^ - - 2 J + 1 

— 2k+l -2”=! £4^^^^^^) +^(2:7=1 ^5 

cos(2; + l)7rct „oo „» sin(2; + l)7rct- 1 ,„oo „* cos(2fc+l)7rs 

^777-2:; = 1 ^6 ^777—) + to 2k+l 

sin(2fc+l)7rs-- ^ __ 1 cosr2 i+ll7rct „„ 


- 2k+l 2”=! £'8 

+ 2c(r, + ^ 0:7=1 £,* - 27=1 Fi*o 

^,'\^oo I-* cos(2k+l)7rs „oo I-* Sin(2fc+l)7rs. 1 ,„oo i-* cos(2; + l)7rct 

)(2:k=i £11 ^777-2:j=i £12 2fe+i ) + 777(^J=1 ^13 2,+i 

- 27-1 £ 1*4 + i- ( 5 ]“ _ 5 ;“ £■* 

2^1-1 14 2j + l ^ 2fc+l 2jk-l 16 2k+l 

))0,(t) + c^(r« + (27=1 £ 1*7 - 27=, El , 


( 2 : 7=1 ^21 


(27=i£i% 2T+; -2:^=1 

^00 r-* sin(2; + l)7rct^ , 1 cos(2k+l)7rs r-* sin(2k+l)7rs 

- Z,; = l C22 7777 ) + — VZ,k=l C23 7777 Z,k=l C24 7777 

))<?q(t)] = Mg(p^„,{ct)^rn{.s) 


On rearranging and simplifying equation (3.59),one obtains 

1 V"" COS(2/ + 1) TTCt 


(1 + VJ(p*(T, + - 


cos(2fc + l)7rs 


2; + l 


Z ^.sin (27 + l)7rct^^'^ , 

._ 27 Ti 


2/C + 1 


sin(2fc + l)7rs^ ^ 1 cos(27 + l)7rct 

2fc + 1 ^ ’’’ 471 ^ A 27 + 1 


Y sin(27 + l)7rct , 

A 27 + 1 
))(?„(t) + 2cm(p*(Ty 


1 V-' 


cos(2fe + l)7rs 


-E' 


1 Y" cos( 27 + l)7rct Y 

7r2 A 27 + 1 ^ 


2; 3 

sin(2fe + l)7rs 
2fe + 1 ^ 

sin(27 + l)7rct 
2 j + 1 
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^^voo n* cos(2k+l)7rs ^oo n* sin(2fc+l)7rs^ , 1 n* cos(2;+l)i 

KU=1 ^11 - lj=l ^12 2k+l ) + = l 

„oo r-* Sin(2j + l)7rct, 1 ,„(jo r.* cos(2fc+l)7rs 

- ^14 —iTTi—) + ak=i ii5 — 

))<?n(t) + (fr ^ ^ i-T 


2; + l 

, sin(2fc+l)7rs 

= l “14 2; + l "" ' 2k+l 2.k=l 16 2fc+l 

))<?n(t) + (f4 + ni<p‘c^(T8 + ;^Q:7 =i £i*7 

/-voo 17* C0S(2k+l)7rs ^oo 17* sin(2k+l)7rs,j 1 ,„oo 17* cos(2; + l)7rct 
I2*k=1^19 2k+l 2*k=lfi20 2k+l >> +47r'''^J = l ^21 2j+l 

VOO r** Sin(2;+l)7ret^ , 1 ^^00 r** cos(2k+l)7rs ^00 r** sin(2k+l)7rs 

- 2*/=l *22 ) + - {U=l *23 2*k=l *24 2k+l 


2J + 1 4.j = l 2, 

J_*yoo £* cos(2; + l) 7 i:ct 

;k+l ^K=i -^u 2k+l ' ' 471 = 1 2;+l 

(2;+l)7ret^ , 1 r** cos(2k+l)7rs ^oo r** sin(2k+i: 

2J + 1 >' + 47r''^fe=l ^23 2k+l 2*k=l*24 2k+l 

)))Q„(t) + UTcp* [(Te +^.(1.7=1 - ^T=l ^2 

sin(2y+l)7rct,^ cos(2fe+l)7rs ^oo sin(2fe+l)7rs,^ ^ 1 17* 

iTTI K2*k=i *3 2*k=i *4 2k+i ^ ^ 

cos(2; + l)7rct „oo i-* sin(2; + l)7rct. 1 ,„oo i-* cos(2k+l)7rs „oo i-* 

^777-2:j=l *6 2;>1 ) + - 


sin(2fc+l)7rs ^ ^ 1 „oo „* cos(2;+l)7rct „oo p* sin(2;+l)7rc 

2k+l + 71 aj = i *9 ^777-lj = l *10 2J>1 

N^Y'oo p* cos(2fc+l)7i:5 ^oo „» sin(2fc+1)715 ,' , 1 „* cos(2j + l)7ict 

Jv2jk=l 11 2k+l ^; = 1 12 2fc+l 471 = 1 2j + l 

- 2 : 7=1 *1*4 ^ (^”=1 - 2“=1 

))0,(t) + C^(T8 + ;^( 2 : 7=1 *1*7 ^ 

^*_^ cos(2k4i)7i5 _ ± ( 5 : 71 ^ £££(|^ 

^00 r-* sin(2; + l)7ict,' , 1 ,,^00 cos( 2 k+l) 7 is 

- z.;=i C22 7777 J + — (,z,k=i C23 ;rr; 

))<?q(t)] = Mgcp^rnictj^mis) 

On further simplifications and n 


^00 17 * sin( 2 ;+l) 7 ict,' 

2*1=1 ^18 21+1 ^ 


-1 - 2J+1 

yoo p* sin( 2 k+l) 7 is 
2k+l ^k=l 24 2fc+l 


_-arrangement, one obtains 

Q.m + 2ot,.-( 7',+- zr., e;. e;, 

cos(2fe+l)7ict yoo „* Sin(2fe+l)7ict , 1 -ypp , cos(2j + l)7ict ypp , 
2k+l 2*k=l 12 2k+l 471''^-' = 1 2; + l '^-'’=1 I"! 

+ (^”=1 ^15 “'^2 kir '' - ^”=1 + 

(fL(l -nr<p*(*a +;^(2:7=i *i* - 2:7=i 

^* cos(2k.l)7icl _ Sinc^^^ ^ ^ 005(2)^ _ 

^00 t7* sin(2;+l)7ict,' 

-^;=1 *6 2J+1 ^ 


E oo sin(2;+l)7ict, ,„oo , cos(2k+l)7ict „(xj „* sin( 2 fc+l) 7 ict,,,, 

1 = 1 21 + 1 ) + <-^fc=l -^fe=l 2 k+l )))) ^ 

c^rn<p*iT,-^^a:T=^ *7 /°-Pi^i)--i _^7i^ *i*8 ""^2i+?"' )g”=i ^1-^ 

... 


2 fc+l 

, 21+7-2^7=1 ^IB—7771^)0 

cos( 2 fc+l) 7 ict „oo I-* sin( 2 fc+l) 7 ict, 1 ,^00 1-* cos(2J+l)nct 

^777 -Zk=l *20 2 fc+l ) + ^ <-^ 1=1 ^31 2141 

""^21>?"' ) + (^“=1 ^33 “'^2ktr'' - 2“=1 ^3*4 
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+^<p* [(T’e + ^0:7=1 ^1* 

cos(2k+l)nct , sin(2fc+l)7rct, 1 

Lk=i ^4 r—; ; + 7z\Li 


1.7=1 

^voo r-* cos(2; + l)7rct 

-a.j=i is 


3s(2k+l)nct yoo sin(2fc+l)7rct, ^ 

2fc+l 2,k=l 4 2fc+i ^ 

Sin(2;+l)7ret^ , r* cos(2fc+l)7rct 


^TTI- 

-I /T I 1 /'\-'oo „» cos(2; + l)7rct „oo i-* sin(2; + l)7rct. ,„oo i-* cos(2k+l)7rct 
2 ciTy + - - Zj=i ilo 2y+i ill 2k+i 

- Sk-l il‘2 + i- (2“ 1 cos(2j + l)7rct _ sin(2; + l)7rct 

2J/C_1 IZ 2fe+l ^ 2/+1 2/+1 ^ 


cos(2y+l)7rct v’oo r** sin(2j + l)7rct,^ j.,* cos(2fe+l)7rct v’oo r** 

^771-^18 2y+i )ak=i ii9 -lfc=l i20 

sin(2fc+l)7rct^ ^ 1 r-* cos(2j + l)7rct ^oo r-* sin(2j + l)7rct^ ^ 1 

^ = 1 ^21 iTTI hj=l i22 ^77^^ J + - 

017=1 i23^^^^g^-2:”=l i24^^^^1^^))Q,(t)] = Mg<p<i>mict)<i>mis) 
Applying modified asymptotic method of Struble, the solution to equation (3.61) becomes 
(?n(t)+>9n"(3n(t) = 0 

Hence, the entire equation becomes 

Qn(t) + ^nQn(t) = M g(p^rn(ct)^m(s) 


where 


= fnn -^(^Ltn(p*(Te +^07=1 E; 


,* cos(2k+l)7rs 


-i:k=ii8 


iin(2k+l)7i 


')) 


2 1 ^x-'oo I-* C0S(2k+l)n:ct „aD r'* sin(2k+l)7rct,,, 

-c^mcp (Tg + -(2:k=i i23-- 1.7=1 i24——))) 

which is the modified frequency representing the frequency of the free system. 

Rewriting equation (3.64), one obtains 

Quit) + dnQnit) = W5<P‘I’m(5)[sin“m(f) + Aj^cosa^t + B^sinha^t + C^cosha^t] 
Following the procedures applied to solve equation (3.51), one obtains 

Qnit) = [(“m + ’9^)(“mSinT9„t - l9„sinamt) - Am'dnicCm + l9^)(C0Sa:mf 

-cosi9„t) - BmioCm - i9^)(amSini9„t - i9„sinhamf) + Em^niam - i9n)(coshcfmt 
-cosi9„t)] 

which on inversion yields 


W(x,y,t) = T.7rn=i I“m=i + >9n)(«mSmT9„t - i9„sinamt) - Am^niai, 

+i9^)(cosamt - cosi9„t) - - i9^)(cfmSini9„t - i9„sinhamt) + C^i^nC^m “ ’^n) 

(cosha^t - cos79„t)](sin^x + Ap^cos^x + Bp^sinh^x + CpmCosh^x) 
(sin^y + Aq^cos^y + B,^sinh^y + C,^cosh^y) 


which is the transverse displacement response to a moving mass of a rectangular plate. 


(3.61) 

(3.62) 

(3.63) 

(3.64) 

(3.65) 

(3.66) 

(3.67) 
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"W. ILLUSTRATIVE EXAMPLES 


I Orthotropic Rectangular Plate Clamped at All Edges 

For an orthotropic plate clamped at all its edges, the boundary conditions are given by 
V(0,y,t) = 0, V(L^,y,t) = 0 
V(x,Q,t) = Q, V(x,Ly,t) = Q 
sviLx.y.t) _ 


Thus,for the normal modes 


SV(0,y,t) _ 

dVix,0,t) _ 
dy 


fpm(O) = 0, 
= 0 , 
afprnCO) _ n 




For simplicity, our initial conditions are of the form 

V{x,y,G) = 


dVix.Ly.t 

dy 


^=0 


= 0 


^bm(Ly) = 

ajymilx') ^ 
dx 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 


Using the boundary conditions in equations (4.1) to (4.4) and the initial conditions given by equation (4.9), it can be shown 

that 


In the same vein, we have 


Apm 


A, 


sinhfpm-sinfpm rosfpm-coshfpm 


cosfpm-coshfpn 

i sinfpn 

x+sinhfpm 

sinhfftm-sinfsm 

^ _ COSfp„ 

i-coshfftm 

cos^tm-cosh^trt 

1 sinfftn 

x+sinhfftm 

sinhf^-sinf^ 

cosfrn- 

■coshf,„ 

cosfm-coshfm 


■sinhfm 



and from equation (4.12), one obtains 

cosfmcoshfp, = 1 

which is termed the frequency equation for the dynamical problem, such that 

= 4.73004, ^2 = 7.85320, = 10.9951 


(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


On using equations (4.12), (4.13), (4.14), (4.15) and (4.16) inequations (3.57) and (3.67), one obtains the displacement response to 
a moving force and a moving mass of clamped orthotropic rectangular plate resting on bi-parametric elastic foundation 
respectively. 


V. DISCUSION OF THE ANALYTICAL SOLUTIONS 

For this undamped system, it is desirable to examine the phenomenon of resonance. From equation (3.57), it is explicitly shown that 


http://www.aipublications.com/iicmp/ 


Page I 84 











Open Access 


International journal of Chemistry, Mathematics and Physics (IJCMP) 
[Vol-4, Issue-4, Jul-Aug, 2020] 
https://dx.doi.ora/10.22161/iicmp.4.4.2 

ISSN: 2456-866X 


the clamped orthotropic rectangular plate resting on constant elastic foundation and traverse by moving distributed force with 
uniform speed reaches a state of resonance whenever 

«m = fnn (5-1) 


while equation (3.67) shows that the same clamped orthotropic rectangular plate resting on constant elastic foundation and 
traverse by moving distributed force with uniform speed reaches a state of resonance when 


where 






for all values of n. 


(5.3) 


Conclusively, for the same natural frequency, the critical speed for the moving mass problem is smaller than that of the moving 
force problem. Thus, resonance is reached earlier in the moving mass system than in the moving force system. 


VI. NUMERICAL CALCULATIONS AND DISCUSSION OF RESULTS 

To illustrate the analysis presented in this work, orthotropic rectangular plate is taken to be of length Ly = 0.923m, breadth = 
0.432m the load velocity c=0.8123 m/s and s = 0.4m.The results are presented on the various graphs below for the simply 
supported boundary conditions. 

6.1 Graphs for Clamped - Clamped End Conditions 

Figures 6.1 and 6.2 display the effect of foundation modulus (Ko) on the deflection profile of clamped-clamped orthotropic 
rectangular plate under the action of load moving at constant velocity in both cases of moving distributed forces and moving 
distributed masses respectively. The graphs show that the response amplitude decreases as the value of foundation modulus (Ko) 
increases. 



Fig. 6.1: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Ko and Traversed by Moving 

Force 
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Fig.6.2: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Ko and Traversed by Moving 

Mass 


Figures 6.3 and 6.4 display the effect of shear modulus (Go) on the deflection profile of clamped-clamped orthotropic rectangular 
plate under the action of load moving at constant velocity in both cases of moving distributed forces and moving distributed masses 
respectively. The graphs show that the response amplitude decreases as the value of shear modulus (Go) increases. 



Fig. 6.3: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Go and Traversed by Moving 

Force 
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Fig.6.4: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Go and Traversed by Moving 

Mass 


Figures 6.5 and 6.6 display the effect of flexural rigidity of the plate along x-axis (Dx) on the deflection profile of 
Clamped-clamped orthotropic rectangular plate under the action of load moving at constant velocity in both cases of moving 
distributed forces and moving distributed masses respectively. The graphs show that the response amplitude decreases as the value 
of flexural rigidity (Dx) increases. 


5E-11 


~ 3E-11 - 

'■A /“» //i 

t 3E-11 

i \ / i / \ / \ 


If-A 

i °' 


I -2E-11 - 

i \ / ' 

1 

S -3E-11 

\.J \} \ji 



-4E-11 

— 

-SE-. 

- Dx=2000 


Fig. 6.5: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Dx and Traversed by Moving 

Force 
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Fig.6.6: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Dx and Traversed by Moving 

Mass 


Figures 6.7 and 6.8 display the effect of flexural rigidity of the plate along y-axis (Dy) on the deflection profile of 
Clamped-clamped orthotropic rectangular plate under the action of load moving at constant velocity in both cases of moving 
distributed forces and moving distributed masses respectively. The graphs show that the response amplitude decreases as the value 
of flexural rigidity (Dy) increases. 



Fig.6.7: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Dy and Traversed by Moving 

Force 
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Fig.6.8: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Dy and Traversed by Moving 

Mass 


Figures 6.9 and 6.10 display the effect of rotatory inertia (Ro) on the deflection profile of Clamped-clamped orthotropic rectangular 
plate under the action of load moving at constant velocity in both cases of moving distributed forces and moving distributed masses 
respectively. The graphs show that the response amplitude decreases as the value of rotatory inertia(Ro) increases. 



Fig. 6.9: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Ro and Traversed by Moving 

Force 
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Fig.6.10: Displacement Profile of Clamped-clamped Orthotropic Rectangular Plate with Varying Ro and Traversed by Moving 

Mass 


VII. CONCLUSION 

In this research work, the problem of dynamic behavior of 
clamped orthotropic rectangular plates resting on 
bi-parametric foundation has been studied. The closed form 
solutions of the fourth order partial differential equations with 
variable and singular coefficients governing the orthotropic 
rectangular plates is obtained for both cases of moving force 
and moving mass using a method of Shadnam et al and 
another solution technique that is based on the separation of 
variables which was used to remove the singularity in the 
governing fourth order partial differential equation and 
thereby reducing it to a sequence of coupled second order 
differential equations. The modified Struble’s asymptotic 
technique and Laplace transformation techniques are then 
employed to obtain the analytical solution to the 
two-dimensional dynamical problem. 

The solutions are then analyzed. The analyses show that, for 
the same natural frequency and the critical speed for the 
moving mass problem is smaller than that of the moving force 
problem. Resonance is reached earlier in the moving mass 
system than in the moving force problem. That is to say the 
moving force solution is not an upper bound for the accurate 
solution of the moving mass problem. The results in plotted 
curves show that as the rotatory inertia correction factor (Ro) 
increases, the amplitudes of plates decrease for both cases of 
moving force and moving mass problems. The flexural 
rigidities along both the x-axis (Dx) and y-axis (Dy) increase, 
the amplitudes of plates decrease for both cases of moving 


force and moving mass problems. As the shear modulus (Go) 
and foundation modulus (Ko) increase, the amplitudes of 
plates decrease for both cases of moving force and moving 
mass problems. 

It is shown further from the results that for fixed values of 
rotatory inertia correction factor, flexural rigidities along both 
x-axis and y-axis, shear modulus and foundation modulus, the 
amplitude for the moving mass problem is greater than that of 
the moving force problem which implies that resonance is 
reached earlier in moving mass problem than in moving force 
problem of clamped orthotropic rectangular plates resting on 
bi-parametric foundation. 
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